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Abstract 

The hierarchical scenario of structure formation, in the frame of the A- 
CDM cosmology, predicts the existence of dark matter (DM) sub-halos 
down to very small scales, of which the minimal size depends on the 
microscopic properties of the DM. In the context of annihilating DM, such 
substructures are expected to enhance the primary cosmic ray (CR) fluxes 
originating from DM annihilation in the Galaxy. This enhancement has 
long been invoked to allow predictions of imprints of DM annihilation on 
the antimatter CR spectra. Taking advantage of the method developed 
by [2], we [T] accurately compute the boost factors for positrons and 
anti-protons, as well as the associated theoretical and statistical errors. 
To this aim, we use a compilation of the latest results of cosmological 
N-body simulations and the theoretical insights found in the literature. 
We find that sub-halos are not likely to significantly boost the exotic 
production of antimatter CRs. 
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Introduction: Like other topics in fundamental physics, annihilating dark 
matter is rather well motivated by some coincidence arguments, because micro- 
scopic new physics is expected to surge at the electroweak energy scale: This 
energy scale gives also the right interaction amplitudes to naturally get the 
correct abundance of WIMPs today to account for the cosmological DM en- 
ergy density, provided that they are thermally produced in the early Universe, 
without matter-antimatter asymmetry (for a recent review on dark matter can- 
didates, see the nice lecture by [3]). The dark matter annihilation properties 
are intimately connected to its possible direct interaction with normal mat- 
ter, and this sketches two complementary approaches in order to hunt for dark 
matter signatures, beside its production in particle colliders. Direct detection 
experiments focus on energy deposits on very sensitive and deep underground 
detectors, while indirect detection aims at observing annihilation traces in the 
form of neutral or charged cosmic rays (for a review on indirect detection, 
see e.g. [4]). Primary antimatter CRs originating from DM annihilation could 
be hunted e.g. in the positron, anti-proton and anti-deuteron spectra at the 
Earth |5j. In particular, the positron spectrum has some features around 10 
GeV that are still hardly understandable from a secondary origin [5] . Some au- 
thors tried to explain this with DM annihilation induced positrons, but most 
of the predictions rely on a rescaling of the primary fluxes by invoking some 
constant boost factors due to DM dumpiness (see e.g. [8j |9]). Indeed, many 
DM substructures are predicted in CDM-type cosmologies, and as the annihi- 
lation rate is proportional to the squared DM density, it naturally increases in 
an inhomogeneous medium, for a given average density. Though the idea of 
enhancement due to sub-halos is far from new [7], it has never been fully ad- 
dressed in the context of charged cosmic rays, contrary to gamma-rays, due to 
the additional difficulty of dealing with CR propagation. Furthermore, because 
dark matter clumps might be numerous and countable (at least theoretically 
or in simulations), they are better described by some phase space distributions. 
The full treatment of this problem has therefore to make use of statistics, and 
should also provide either predictions and the associated statistical uncertain- 
ties (beside the theoretical ones). As a matter of fact, from an observational 
point of view, any signature can be said unambiguous if, at least, the theo- 
retical predictions do not suffer from large uncertainties, unless the expected 
signature is very specific and clear. [2] developed a detailed method in order to 
translate the dumpiness phase space in terms of probability for the primary sig- 
nal/boost. These authors focused on positrons, and took a very simple model 
in which all sub-halos had the same properties and were spatially distributed 
according to the host DM halo profile. They showed that the boost factor and 
the statistical uncertainties actually non trivially depends on the detected CR 
energy. We have used their method, which is quite general and applies to any 
CR species, but with more sophisticated and more precise inputs for the dark 
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matter distribution. This proceeding summarises the work presented in [I], to 
which we refer the reader for a more exhaustive bibliography, and in which we 
accurately compute the boost factors resulting from dark matter dumpiness 
for positrons (e + 's) and anti-protons (p's). 

Dark Matter distribution in the Galaxy: Any prediction for indirect 
detection of DM strongly depends on the DM distribution that is used, for the 
smooth halo as well as for sub-halos (we refer to the smooth dark matter com- 
ponent as the host DM halo of the Milky Way — MW). The density profile of 
dark objects is the main quantity to anticipate how the DM annihilation rate 
may spatially evolve in the Galaxy. Though axysymmetry characterises most 
of the results of cosmological N-body simulations, we adopt here a spherical 
symmetry to describe the Galactic host halo and the sub-halos. A structure of 
any scale may thus have its density profile p(r) depicted by the generic formula 
proposed by [TU]: p(r) = p s / [(r/r s ) 7 [1 + {r/r s ) a ]^ l3 ~ 1 ^ al \ , where p s is twice 
the density at the scale radius r s , at which the logarithmic slope changes from 
7 to ((3— ~f)/a. The well-known NFW [TT] and Moore [T2] profiles are recovered 
with {a,/3, 7} — {1,3,1} and {1.5,3,1.5} respectively. An NFW behaves like 
r _1 while a Moore scales like r -15 in the central region of the structure. The 
scale radius is connected to the concentration parameter that we will discuss 
further. The DM density may saturate at the very centre of a structure because 
of annihilation, and a cut-off radius is usually set by equating the gravitational 
infall rate with the WIMP annihilation rate, of which the typical size is ~ 10~ 6 
pc for usual WIMPs. For the smooth component, we choose r s = 20 kpc, 
and we normalise the density profile p sm with respect to the local DM density 
in the solar neighbourhood p sm (R,Q — 8 kpc) — p& = 0.3 GeV/cm 3 . The 
inner shape of the smooth component has no significant impact on the DM 
annihilation induced CRs because of propagation effects that strongly dilute 
what originates from the central regions of the Galaxy. Fluxes at the Earth 
are much more sensitive to modifications of local distribution of DM. Finally, 
we set the Galactic virial radius to R^i r = 280 kpc, so that the total halo mass 
inside R% it is M^ h . = 1.12 x 10 12 M Q . Sub-halos are found to b c more and more 
numerous as the resolution of N-body simulations gets thinner and thinner. 
For instance [T3] found ~ 10 15 substructures of Earth-mass and of solar system 
size, in a galaxy-size box at early stages of structure formation (z = 26). Such 
small masses are consequent of very small free streaming scales for the dark 
matter, and are rather generic for WIMPs. However, one can wonder whether 
or not such small structures may survive tidal effects and encounter events 
that characterise their evolution in the host halo. While this issue is still de- 
bated, we have scrutinised the effect of considering different minimal masses — 
10~ 6 , 1 and 10 6 M — in the original paper. A complete sub-halo phase space 
may be portrayed with the following normalised probability density functions 
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(pdf), which define the number density of sub-halos per mass unit in the MW: 
%${2,M) = ^% = N x x , where we assume that the 

mass function has no spatial dependence. The mass function is usually found 
to be a power law ; we therefore define: dVl ^ M M ^ = KmM -0 " 1 , where Km 
normalises the pdf to unity in the mass range used for sub-halos, so that it 
depends on M m ; n , Af max and a m . For simplicity, we have considered the same 
power law all over the mass range, which is not that obvious but motivated 
by theoretical arguments [E] and simulation results [T3]. Regarding the spa- 
tial pdf, we have retained two cases. The first one is that sub-halos spatially 
stick to the smooth host DM halo density profile, which may be an appropriate 
description for very light clumps, as they may behave like test particles. In 
that case, dVv(x)/dV — /9 sm (x)/M^ ir , which is normalised inside R^ ir . Nev- 
ertheless, according to the current results of N-body simulations for which the 
resolved sub-halo mass is > 10 6 M© (c/. e.g. the Via Lactea [14 ), the spa- 
tial distribution of substructures is mostly anti-biased compared to the smooth 
component. This is characterised by a spherical isothermal distribution with a 
core radius r c roughly equal to the scale radius r s of the MW. This is the 2 nd 
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case: = R v x 
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, where Ky normalises the distribution to 



unity within the virial radius R\ Y of the MW. The total number of sub-halos 
No depends on a m and M m ; n , but there are constraints from current N-body 
simulations. Indeed, there is a rather good agreement about the counting of 
well resolved sub-halos inside a MW-type host: ~ 100 objects are found in the 

mass range 10 8 -10 10 M Q . Therefore, we ask for N x J^m^ dV M /dM = 100. 
The inner profile p c \ of any sub-halo is also described with the NFW or the 
Moore models, the latter defining a maximally optimised scenario. Besides, 
we need additional prescriptions in order to set the scale radius r s and the 
associated scale density p s for sub-halos. The standard method is to define 
some reference quantities and to connect them with the physical scale variables 
by some functions usually fitted on the N-body simulation results. Those ref- 
erence quantities are the virial mass and radius, which are related through: 
Af v i r = ^ri?vir 3 x (A v ; r (z = 0)Qm(z = 0)p C rit(z = 0)). Such a definition is not 
physical in the sense that any mass is related to a radius enclosing A v j r times 
the background matter density today (we set A v i r =340). Indeed, sub-halos do 
not evolve in a constant background, but inside host halos where the density 
is not the background density. Moreover, the previous equation does not carry 
the formation history of the structure, e.g. the fact that lighter clumps should 
be denser (because formed in a denser universe). This is actually encoded in 
the concentration parameter c V i r , which connects the virial radius, as defined 
above from the mass, to the more physical scale radius: c V i r = -p 2 ^ , where r_2 
is the radius at which d/dr (r 2 p c i(r)) |r=r_ 2 = ( r -2 — r s for an NFW profile). 
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There are different concentration models in the literature, and we have taken 
two extreme cases that encompass them (i) the — maximal — Bullock et al [H] 
model (B01) and (ii) the — minimal — Eke et al [H] model (ENS01). Once the 
sub-halo mass and the mass-concentration relation are known, all the clump 
properties are fixed, and one can compute the corresponding annihilation rate. 
As it is proportional to the squared DM sub-halo density, it is useful to define an 



effective annihilation volume for any sub-halo as follows: £ = J , d 3 x 



where pq , the local dark matter density, allows for a normalisation to the local 
annihilation rate. Such an effective volume is that of the DM within the clump 
would have if it were diluted down to the density p Q . It is convenient to nor- 
malise to local quantities because cosmic ray propagation favours a local origin 
for the exotic contribution. Therefore, £ = f(M), and the mass pdf oIVm/oIM 
translates very simply to an annihilation rate pdf. Changing the inner profile 
results in a shift by a constant factor, e.g. £ m00 rc — 10 x £ n fw, whereas changing 
the concentration model is not as straightforward. Nevertheless, we find very 
roughly that £boi — 1 — 10 x £enso1j depending on the sub-halo mass. 

Cosmic ray propagation: The CR propagation modelling is a key ingredi- 
ent for this kind of studies. Here, we adopt a slab diffusion zone, featured by 
its radial extension that we fix to i? s iab = 30 kpc, and by its half-thickness L, 
3 kpc here. CRs are either confined within the slab or escape forever, which 
is merely fulfilled by imposing Dirichlet boundary conditions to the diffusion 
equation (the CR number density vanishes on the borders). Regarding the 
transport processes, the spatial independent diffusion coefficient is given by 
K(E) = f3K lZ s (where 1Z = pc/Ze is the rigidity) and a constant convective 
wind y conv is directed outwards along the vertical axis. Such a configuration is 
quantified with the medium set of parameters provided by [18 : Kq = 0.0112 
kpc 2 /Myr, S — 0.7 and V conv = 12 km/s. One can easily write and solve 
the diffusion equations for both e +/ s and p's for this kind of geometry. For 
e + 's, the main processes that come into play are the energy losses (mainly in- 
verse Compton diffusion off CMB or IR photons, and synchrotron radiation), 
and the diffusion on the magnetic turbulences. Disregarding the convection 
process, which is much less efficient than energy losses, one can express the 

/ / e S ~ 1 —e S ~ 1 \ \ 1 / 2 

typical propagation length for e + 's as: Ad = I 2KqTe I — jrj 2 — ) ) 7 where 
e = E/{Eq = 1 GeV}. Assuming an infinite 3D spherical diffusion zone (cor- 
rect while Ad ^$ L), the e + propagator is proportional to a Gaussian function 
of the source distance with a = Ad (< few kpc): Sources located farther than 
Xd will almost not contribute to the flux at the Earth. Ad being a decreasing 
function of the detected energy, the effective volume in which e + 's propagate 
increases as they loose energy. The propagation of p's must include spallation 
processes and wind convection that occur in the thin Galactic disc, so we can 
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not use a simple spherical symmetry to derive a global expression. Moreover, 
the energy losses are negligible for this species, which modifies significantly the 
picture that we had for e + 's. Nevertheless, it is useful to write, as for e +/ s, the 
typical propagation length: Ad = , where convection is assumed to dom- 

inate over spallation in average, which is correct unless at sub-GeV energies. 
This is quite different from e +/ s because this length is an increasing function 
of energy. The picture is therefore reversed, and the propagation volume is 
much larger at higher energy for p's (Aj reaches the size of the diffusion slab 
at energies of order 10-100 GeV). We now define a convenient Green function 
for any CR species, that encodes the injected spectrum induced by dark mat- 
ter annihilation: G(E,x Q <- x) = dE s G(E,x Q <- E Sl x) x dN f^ s) , 
where dNcn/dEs is the injected spectrum at source (E = E$ for p's). 



Exotic fluxes and associated boost factors: The Galactic host halo is 
described by a smooth dark matter distribution p sm , so that the corresponding 

primary CR flux reads: 4> sm {E) = j^S J hSilo d 3 x G(E,x Q <- x) where 

v is the CR velocity and S = S (a ann v) pq 2 / (2m x 2 ) encodes the main WIMP 
propertied Sub-halos can be considered as point-like sources, and the flux due 



wo 



to the i object is merely: 4> C \^{E) = ^ x S X & x Q(E, x Q <— x{). Then, 
have to sum over the whole population. We can derive a statistical prediction 
by integrating over the sub-halo phase space. The overall clump contribution 
is thus: cj) d , tot {E) = ^o/dM((M)^x Jd 3 xG(E,x Q «- x)^^ = 
-Wo(<£ci) = j^SN n ({;)(g} where the 1 st equality takes directly the — normalised 
— pdfs into account (this limit would be reached e.g. for an infinite number of 
MC realisations). Those pdfs characterise £ and Q, the former being a function 
of the sub-halo mass and the latter being spatially weighted with dVv /dV. The 
last equality gives the same quantities in terms of statistical mean values. This 
assumes no correlations between the considered variables. The calculation of 
the variance of CR fluxes originating from sub-halos <J c i,tot is straightforward: 



T d.tot _ 1 / °£ , °g , °S g g 



where u x is the variance corresponding 



to any variable x. The boost factors for e + 's and p's are the ratios of the CR 
fluxes originating from a clumpy halo to those calculated for the host smooth 
halo alone. Normalising the whole dark matter average density profile at the 
Earth, the effective boost is: B e g(E) = (1 — / ) 2 + </> c i,tot/</>sm , which depends 
on energy. / Q is the average local DM density in form of clumps (the smooth 
component density is (1 — /o)p sm when sub-halos are added). Estimates of 
the mean fluxes for e + 's and p's and associated statistical variances have been 



2 <5 = 1/2 for Majorana particles, 1 otherwise. 
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performed by using the semi-analytical method proposed by [5] ■ In order to ex- 
haustively characterise the effect of all the relevant variables in the problem, we 
have varied the minimal mass of sub-halos M m j n , the logarithmic slope of the 
mass function a m , the concentration model, the inner profiles of sub-halos, their 
spatial distribution in the MW ; we have also varied the CR propagation model. 
This provides a way to anticipate predictions for any dumpiness modelling, and 
is also useful to sketch some theoretical uncertainty contours. The full results 
are available in the original paper pQ. In this proceeding, we quickly discuss 
three cases: minimal, reference and maximal models, considering a medium 
set of CR propagation parameters. By defining a model with {M min , a m , 
inner profile, concentration model, spatial distribution}, the minimal reads 
{10 6 M Q ,1.8,NFW,ENS01,cored }, the reference {10- 6 M Q ,1.9,NFW,B01,cored 
} and the maximal {10 _6 M Q , 2.0, Moore, B01, smooth-tracking }. The most op- 
timistic — maximal — case, which is also the most unlikely, yields boost fac- 
tors ~ 20, which depend on the CR energy, while the other configurations do 
not permit substantial enhancement, the effective boost factors B c g remain- 
ing close to unity for both e + 's and p's. Taking back NFW profiles in the 
maximal setup, we would obtain ~ 20/10 = 2. The variance affecting those 
predictions decreases with the local average number density of sub-halos. In- 
deed, the statistical variance depends on the number of objects located inside 
a volume bounded by the CR propagation length and contributing to the flux 
at the Earth. As the propagation length is a function of the CR energy which 
depends on the species — increasing (respectively decreasing) with energy for 
p's (e + 's) — one can understand why the picture can be different for different 
CR species. As a very general statement, the variance on primary fluxes is 
smaller at lower energy for e + 's, or at higher energy for p's. Nevertheless, even 
with a large variance for sub-halo fluxes, the statistical error will still be small 
for boost factors as soon as the smooth contribution will dominate (i.e. errors 
are small for small boost factors). 

Conclusion: In pQ , we have tried to exhaustively tackle the problem of DM 
dumpiness effects in the frame of indirect detection of annihilating DM with 
antimatter cosmic rays. We have taken the whole phase space of sub-halos 
into account, using a compilation of results found in the literature, coming 
from N-body simulations as well as from analytical models. We have sketched 
the theoretical uncertainties affecting the useful dark matter parameters, and 
defined rather wide associated ranges to account for them. For the sake of 
completeness, we have also considered different sets of CR propagation param- 
eters, which are still degenerate and not enough constrained by the existing 
secondary /primary CR measurements. We refer the reader to the original pa- 
per for very detailed results. In this proceeding, we have only illustrated the 
cases of extreme configurations, which can provide boost factors < 20 for both 
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e+'s and p's. Nevertheless, our analysis strongly disfavours large and even 
mildly boost factors for most parts of the large parameter space used to de- 
scribe sub-halos. This is almost independent of CR propagation uncertainties. 
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